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Abstract. 

This study describes the surface deformation in lower Jehlum and eastern 
Potwar using remotely sensed Shutter Radar Topographic Mission (SRTM) Digital 
Elevation Model (DEM) to identify regional uplift rates. Mapping of active surface 
deformations help to predict seismic hazards. DEMs are the fundamental input factors 
that identify the surface geometry and its belongings. These belongings include stream 
Offset, elevation, and slope breaks within a contributing area. This study provides an 
evidence of tectonic activity and its impact on regional drainage network using SRTM 
DEM. Various indices including concavity and steepness were computed using power law 
in steady state conditions. We prepared the drainage network map of the study site 
showing uplift rates in mm/year. The deflection in stream network proves the 
existence of active fault in this region which controls the local drainage network. The 
results prove the relative uplift along Main Boundary Thurst (MBT) and the impact of 
active tectonic on evolving young organs. 
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Introduction. 

The controls on tectonic structures in a hilly area are technically sound 
regarding horizontal translations, rock uplift, bed rock erosion, climatic variations and the 
development of incised structures [1]. River pattern demonstrates multiple forces in 
action e.g. river geometry is useful to identify the spatial patterns of uplift in variouscases 
[2]. Longitudinal river profile can retort to climate by erosional and hydrological 
situations and also affected by orographic precipitations [3,4]. 

Lithology is a fundamental element which defines the concavity of the river[5]. 
Lithology of a river describes the rate of erosion and incision. This study describes the 
neo-tectonics in the investigation site e.g. faults in eastern Potwar and the stream 
interactions. It is based on remotely sensed analysis of tectonic activity and its impactson 
regional surface deformations. Digital Elevation Model (DEM) is used here to define 
the surface geometry in details. 

Tectonic movement along fault is characterized by various features including 
stream offset, elevation, moraines, landslides, slope breaks, fluvial terraces and the 
drainage zone. We performed the stream profile analysis to compute surface 
deformation in the study site. Stream analysis is important to draw stream network and 
their characteristics in the study site. This analysis is useful to study tectonic activity 
indirectly. 

This study aim at computing steepness, and concavity which is helpful to 
compute the uplift rates. It also aims at investigating the impacts of Muzaffarabad 
earthquake 2005 in Pakistan. 

Study Site and Tectonic Settings 

The study site lies between 73-74E longitude and 33-34N latitude. A part of 
Main Boundary Thurst (MBT) is lying underneath the study site in the northwest 
direction. Various structural styles are located underneath the study site e.g. Riasi and 

Koth Thurst is in SE-NW oriented and the Jehlum fault is also in the same orientation. 
Another fault “Rawat” is oriented in SW-NE direction. The lineaments and out crops are 
tight angle to the direction of transportation at a rate of 17-19 mm/year [6]. Thestudy 
site is mapped in the figure 1: 
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Figure 1. Tectonic map of the Hindu Kush-Himalaya-Pamirs-Karakoram showing 
reported and newly confirmed faults with inset showing the study area. 
Fieure 2 is showing that topography is more heave in NE part in comparison to 
southern patt. 
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Figure 2. Location of study area with earthquakes locations in red and thick black lines 
represent thrust faulting with teeth symbols showing direction of thrust movements. 
Data and Methods 
SRTM Digital elevation model (DEM) is used in this study to extractdrainage 
network and to compute steepness, slope, and concavity. 
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Figure 3. Flow chart showing various algorithms implemented during SRTM DEM 
processing. 


DEM based Drainage Extractions 

D8 algorithm is commonly and mostly used algorithm to compute flow 
direction, accumulation, stream definition, stream segmentations, catchment grid 
delineation, catchment polygon processing and the drainage line processing. A 
depression less DEM 1s required as basic input to compute the flow direction. We used 
the built-in utility in Arc GIS 10.1. “Fill” to generate a depression DEM. This 
algorithm check depressions pixel by pixel and fill them by interpolating nearby pixel. To 
compute “flow direction”, the elevation of central cell is compared with its eight 
surrounding cells and the direction leads to the cell having lowest elevation. In the 
proceeding step, flow is accumulated to define streams. The density of streams in 
drainage network is defined depending upon user’s requirement. If user need high 
stream density, a small number of cells is gtven in input and vice versa. In the next 
step, streams are segmented by connecting points to define the stream order. The 
stream order is useful to compute catchments areas and to draw drainage network in 
vector format. The river was demarcated in the drainage network using river tool 2.4. 
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Figure 4. The diagram showing mechanism of River profile development, modified 
fromSnyder et al (2000). 
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Concavity and Steepness Indices 

In a diverse topography, the DEM shows scaling, where the stream’s slope is 
expressed by a function of Stream power law (SPL) [7,8] which is defined with formulaas 
below 


s=ksA-6@ 

Where A is the drainage area of upstream, 8 and Ks are concavity and 
steepness indices respectively, and “s” is the slope of particular channel. Index of 
concavity (LOC) is represented by@. 

Results and Discussion 

A consistency was not observed in stream river channels due to diverse neo- 
tectonic activities in NE and NNE parts of the study site. The pattern of spatially 
distributed stream channels represents local deviations in thurst belts. The information of 
stream profiles had detailed information of neo-tectonics signals. 

The inconsistency in drainage profiles demonstrates the variation in uplift rates in thestudy 
site. As shown in figure 5. 
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Figure 5. River profile analysis for the Jhelum and Poonch Rivers, three knickpoints are 
identifiable. 

In a watercourse, a change within two index of steepness (LOS) is generally 
connected by a higher or lower index of concavity (IOC) value. However, [OC is 
comparable to early and ultimate profiles of drainage network [9]. The lower value of 
IOC describes the variability in rock strength. The degree of roughness in knick point, 
describes the tectonic events, and lithological disparity or stream capture events. 

A sharp knick zone is an indication of most recent tectonic activity [10]. The 
stream no: 1 in figure 5 of Jehlum river profile shows the IOS value=532.54 and 
IOC=0.48. These values define that area corresponding to stream 1 is less eroded and 1s 
in uplifting process. KP1 describes Kotl Thurst (KT), KP2 represents Quaternary 
Fault (QF2), and KP3 indicates the interaction of Ponch river with Quaternary Fault 
(QF3). KP1 and KP2 is related to intersection of Jehlum Fault (JF) with Jehlum River 
KP3 is under the influence of QF3 in north of Mirpur which proves that Jehlum River is 
under the influence of JF. 
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The sudden variations in geomorphic expressions of stream profiles describes 
variations in horizontal or vertical tectonic movements in lithological contrasts, high 
values of Hack Gradient Index (HGI) were observed in NE and NW part of the study 
site as shown in figure 6, because IOS values are in direct relation to Relative Uplift 
Rates (RURs). 
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Figure 6. The interpolated Hack SL-gradient map for the study area. 
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We extracted 109 streams using SRTM DEM and computed IOS and IOC 
values as shown in figure 7 and figure 8 and applied geomorphic indices using stream 
power law (SPL) to all streams. 
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Figure 7. The interpolated [OC map of the study area. 


Jan 2020 | Vol 2| Issue 1 Page | 26 


OPEN ACCESS . - ‘ ; 
B International Journal of Innovations in Science & Technology 


73°O'E 73°30'E 74°0'E 


34°0'N 


33°30'N 


ss __.-+ 
73°30°E 74°0'E 


73°0'E 

Steepness index coding -—4 Faults @ Earthquakes m@ Cities . 

~ High : 1215 0 75 15 30 A 
Km 


==3 Low: 52 


Figure 8. The interpolated [OS map of the study area 
A telative uplift note map was generated for the study site in Mat lab 
environment. The RUR map show a diversity in uplift rates in various parts of the 
study site. 
RURs is between 0.53-1.21 mm/year in NNE of the study site while RURs 
were observed between 0.05-0.51 mm/year as shown in figure 9. These rates were 
between 0.06-0.36 mm/year in the southern section. 
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Figure 9. Map showing RURs of the study area. 
These rates clearly describe that the NNE region of study site is more uplifted in 
comparison to SW patt. 
The highest value of RURs was observed near to MBT between JF and KT. 
This research indicates that the flow direction and the orientations of streams 1s 
influenced by local and regional faults. 
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Most of Jehlum River is controlled by JF and some of its path is under 
lithological control which again run under Manela fault upto Manela reseruior. Local and 
lineaments ate responsible for growth and development of spatial drainage 
distribution at local scale. 

Conclusion 

The automated drainage network extracted using SRTM DEM is significant to 
investigate SPL which provides information regarding active surface deformation and the 
topographic development in a study site. Geomorphic indices of Poonch and Jehlum 
River are important to examine the surface activity across a study area. This study 
reveals that geomorphometric features are effective indicator to study neo- tectonics 
in young topographies with variable uplift rates. 
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